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ABSTRACT 

We present new MIDI interferometric and VISIR spectroscopic observations of 
HD113766 and IID172555. Additionally we present VISIR 11/im and 18/xm imaging 
observations of IID113766. These sources represent the youngest (16Myr and 12Myr 
old respectively) debris disc hosts with emission on ^lOAU scales. We find that the 
disc of IID113766 is partially resolved on baselines of 42-102m, with variations in reso- 
lution with baseline length consistent with a Gaussian model for the disc with FWHM 
of 1.2-1.6AU (9-12mas). This is consistent with the VISIR observations which place 
an upper limit of 0'.'14 (17AU) on the emission, with no evidence for extended emission 
at larger distances. For HD172555 the MIDI observations are consistent with complete 
resolution of the disc emission on all baselines of lengths 56-93m, putting the dust 
at a distance of >1AU (>35mas). When combined with limits from TReCS imaging 
the dust at ~10/J,m is constrained to lie somewhere in the region 1-8AU. Observa- 
tions at '-^IS/im reveal extended disc emission which could originate from the outer 
edge of a broad disc, the inner parts of which are also detected but not resolved at 
10/zm, or from a spatially distinct component. These observations provide the most 
accurate direct measurements of the location of dust at 1-8AU that might originate 
from the collisions expected during terrestrial planet formation. These observations 
provide valuable constraints for models of the composition of discs at this epoch and 
provide a foundation for future studies to examine in more detail the morphology of 
debris discs. 
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1 INTRODUCTION 

Since first being identified using data from the IRAS satel- 
Ute, the debris disc phenomenon has been the subject of in- 
tense study. Spitzer surveys have confirmed that such discs, 
thought to be the debris material left over at the end of 
the planet formation process , are present around ~15% of 
nearby stars (see e.g. IWv8it3l2008l and references therein). 
The spectral energy distribution (SED) of the dust emis- 
sion in most debris discs peaks at >60/im, implying that 
the dust is cool (<80K), and resides in Kuiper belt-like 
regions (^lOAU) in the systems. However, some systems 
have hot dust on scales <C10AU. The current tally of known 
hot debris discs is ^20 across spectral types A-M (13 cat- 
alogued in lWyatt et al.ll200 7a b. and 7 recently discovered 
with AKARI - AFunwaTa et ah 2009 '). This tally includes sys- 
tems with multiple components where the hottest compo- 
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nent lie s at <10AU such as 77 Tel l|Smith et al.ll2009al 'l and 
13 Leo (|Stock et al.ll2010l ). For some multiple dust compo- 
nent systems (e.g. 77 Corvi. lSmith et al] |2009b' and references 
therein) it is suggested that the hot dust component may be 
fed by a parent planetesimal belt coincident with the cooler 
dust belt. For systems without a known cold dust compo- 
nent, alternative models for the origin of the hot dust must 
be sought. 

Amongst these sou rces, HP 113766 ( ^16Myr old) and 
HD172555 (~12Myr old. lZuckermanll200li 'l are the youngest. 
They also have so me of the brightest levels of excess 
IWvatt et al.l l2007al lbl: see Table [1] for further source de- 
tails). The favoured interpretation for the emission observed 
around these sources is t hat we are witnessing ong oing ter- 
restrial planet formation l|Kenvon fc Bromley||200j 'l. since a 
detailed analysis of the Spitzer IRS spectra of both sources 
indicates that the dust composition is similar to that ex- 
pected fr om the catastrophic disr uption of terrestrial planet 
embryos (|Lisse et al. I l2008l . I2OO9I ). However, alternatives to 
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Table 1. Characteristics of science and calibration targets 



Science targets 



Source 


Spectral type Age 


RA 


Dec 


at 10/^m i^jisc a-t lOfim 


Stellar angular size 


Predicted disc size 


HD 


Gyr 






mJy nijy 


mas 


mas 


113766 


F3V 16" 


13 06 35.8 


-46 02 02.01 


94 2359 


0.048±0.003 


13, 31-69= 


172555 


A5V 12' 


18 45 26.9 


-64 52 16.53 


721 973 


0.27±0.01 


198'* 


Standard stars 


Source 


Spectral type 


RA 


Dec 


Fi, at 10/im 


Angular size 


Instrument 


HD 








mJy 


mas 




111915 


K3.5III 


12 53 06.91 


-48 56 35.93 


13797 




VISIR 


110253 


K2III 


12 41 09.67 


-44 06 04.27 


1065 


0.87±0.02 


MIDI 


112213 


MOIII 


12 55 19.43 


-42 54 56.50 


14542 


3.16±0.02 


MIDI 


116870 


K5III 


13 26 43.17 


-12 42 27.60 


10416 


2.58±0.01 


MIDI 


152186 


Kim 


16 55 34.43 


-60 40 38.77 


785 


1.00±0.02 


MIDI 


156277 


K2III 


17 21 59.48 


-67 46 14.30 


7004 


2.00±0.01 


VISIR, MIDI 


169767 


G9III 


18 28 49.86 


-49 04 14.10 


9006 


2.15±0.01 


MIDI 


171212 


Kim 


18 36 41.43 


-56 13 37.12 


754 


1.54±0.02 


MIDI 


171759 


KOIII 


18 43 02.14 


-71 25 21.20 


12640 


2.68±0.01 


TReCS, MIDI 



° Age taken from Sco-Cen association members hip. Age f r om m embership of /3 Pic moving group. From fit to IRS spectra by 
iLisse et al.l | |2008|) . From fit to IRS spectra by iLisse et al.l | |2009| ) . Estimated values of arise from fitting a Kurucz model 
photosphere of appropriate spectral type to the 2MASS K band measured photometry. For the science targets -Fdisc is determined from 
the Spitzer IRS spectra of the target after subtraction of the photospheric model. For the standard targets used in MIDI observations 
angular size was taken from the CalVin tool available at http://virwvir.eso.org/intruments/midi/tools where available (HDs 112213, 
116870, 156277, 169767, and 171759). For the remaining standard targets and the stellar components of the science targets the angular 
size was estimated by assuming that the stars have a diameter typical for their spectral type (taken from Coxll2000i ) and using the 
Hipparcos parallax to determine their distance. Standard stars are used as calibrators for the instruments listed (see text for details). 

present VISIR imaging of HD 113766 together with VISIR 
spectroscopy and MIDI observations of both HD 113766 and 
HD172555. We also present a re-analysis of archival TReCS 
imaging of HD 172555. These observations are compared 
with models for the distribution of the dust and constraints 
placed on the location of the emitting material. 



2 OBSERVATIONS WITH 8M INSTRUMENTS 

2.1 VISIR imaging and photometry of HD113766 

High resolution imaging with 8m-class telescopes can reveal 
debris disc structure on >0'.'5 scales which corresponds t o 
tens of AU for nearby stars (see e.g. ISmith et al.l l2009al l. 
We used VISIR on the VLT to search for emission around 
HDl 13766 on such scales. Observations were performed in 
filters SiC (Ac = 11.85/im, AA = 2.34/im, hereafter referred 
to as N band) and Q2 (Ac = 18.72^im, AA = O.SS^im, here- 
after referred to as Q band) under observing program 079. C- 
0259(A) in April 2007. Observations were performed in a 
perpendicular chop-nod pattern. The science observations 
were calibrated using observations of HD11191 5 taken from 
the Cohen catalogue of mid-infrared standards ( Cohe n et al.l 
1999). Standard star observations were performed immedi- 
ately before and after the science observations to allow mea- 
surements of variations in photometry and the PSF, cru- 
cial in the search for extended emission. A summary of the 
observations is given in Table [21 Data reduction was per- 
fo rmed with cu s tom r outines, the details of which are given 
in ISmith et al.l (|2008h . Data reduction involved determina- 



the planet formation origin model for such hot emission do 
exist, includ ing the sca ttering of comets from tens of AU in 
the syste m jComes et a l. 2005), the sublimation of one su- 
percomet (jBeichman et al..,2005'). a rece nt collision between 
two massive asteroids (ISonget al.ll2005l). a disc of planetes- 
imals on highly eccentric orbits ( Wvatt et al-lfioiol l. or that 
it is in fact a steady-state phenomenon . This has recently 
been suggested for the HD69830 system (|Henell201ll '). which 
had previ ously been identified as a host of transient debris 
emission, (jWvatt et al.ll2007al \ 

A correct interpretation of the hot dust depends crit- 
ically on its radial location. We expect different dust dis- 
tributions from the different theories for the origin of the 
dust, e.g., a population of comets scattering inwards would 
be expected to be observed at multiple distances from the 
star (the parent belt location and the dust sublimation ra- 
dius) whereas dust from terrestrial planet formation would 
be expected to be confined to a narrow ring. Modelling the 
SED provides poor constraints on the radial distribution of 
the dust, as one may either: (i) underestimate the size of 
the dusty region, because the emitting grains are small and 
hotter than blackbody (observed discs ca n be a factor of 
3 lar ger than predicted due to this effect, ISchneider et al.l 
I2OO6I ): or (ii) overestimate the size of the disc, because the 
dust is in an extended distribution (for example the pre- 
dicted size of the C, Lep disc was more than double the 
observed size as multiple disc components over a range of 
distances fro m the star had been f i t by a single disc tem - 
perature, see iMoerchen et"al] l2007l : ISmith fc WTOt^l2010l ). 
The most direct way to resolve these ambiguities is with 
very high spatial resolution observations. In this paper we 
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Figure 1. The final co-added images for HD111915 (the standard star) and HD113766 observed with VISIR (two leftmost columns). 
The contours are at 10, 25, 50 and 75% of the peak and show that there is no evidence for extended emission beyond the size of the 
PSF. The known binary companion HD113766B is difficult to see in the final images of HD1136766A. In the residual image, created 
by subtracting the final standard star image from the HD113766A image after scaling to the peak, there is no evidence for significant 
residual emission which might be a resolved disc. The signal seen in the N band residual image is from the binary companion. In these 
images black pixels are <lcr per pixel. For the residual image the maximum pixel value (white) is 3a- per pixel. The HD113766B image 
for the Q band contains no signal as there is no detection of the binary companion at this wavelength. 



tion of a gain map using the mean values of each frame to 
determine pixel responsivity (masking off pixels on which 
source emission could fall, equivalent to a sky flat). In ad- 
dition a dc-offset was determined by calculating the mean 
pixel values in columns and rows (excluding pixels on which 
source emission was detected) and this was subtracted from 
the final image to ensure a flat background. Pixels show- 
ing high or low gain, or those that showed great variation 
throughout the observation, were masked off. Finally, the 
four images (two negative, two positive) resulting from the 
perpendicular chop-nod pattern were co-added to give a fi- 
nal image for HD113766 and observations of the standard. 
The center of each image was determined through fitting 
with a two-dimensional Gaussian, with the center of each 
image aligned in the co-addition step. Aperture photome- 
try centered on the Gaussian peaks of the final images was 
performed using l"radius apertures (noise levels were de- 
termined from an annulus with inner radius 2" and outer 
radius 4"). The calibrated flux found in the final images 
was 1673±42mJy for HD113766A at N and 1895±34mJy 
at Q. These uncertainties include calibration errors of 2% 
and 6% respectively determined from variation in calibra- 
tion factors between the two standard star observations in 
each filter. These fluxes are compatible within the errors 
with the IRS photometry, which taken over the filters used 
here are 1599mJy at N and 1867mJy at Q. 

The final VISIR images of the PSF reference (standard 
star image) and HD 113766 are shown in Figure [T] (first and 
second columns respectively; top row is N band images, Q 
band images are shown in the bottom row). The final im- 
ages are shown with contours at 10, 25, 50 and 75% of the 
peak. There is no evidence that there is extended emission 
around HD 113766 in either the N or Q band images from 
the contour lines (as compared to those of the standard star; 
low level ellipticity in the N band standard star image is be- 



Table 2. Log of the VISIR imaging observations of HD113766 
and standard star HD 111915. 



OB ID 


Filter 


Int. time (s) 


Object type 


Star 


265452 


Q2 


400 


Cal 


HD111915 


265453 


Q2 


1800 


Sci 


HD 113766 


265450 


Q2 


400 


Cal 


HD111915 


265455 


SiC 


105 


Cal 


HD111915 


265457 


SiC 


900 


Sci 


HD113766 


265456 


SiC 


105 


Cal 


HD111915 



low a 3(T significance level). To search for low-level extended 
emission the standard star image was scaled to the peak of 
the HD 113766 image (in each band) and subtracted from 
the science image. The result is shown in the third column 
of Figure [T] (labelled residual) . The source appearing to the 
North- West in the N band is a known binary companion dis- 
cussed in the following paragraph. Excluding this source the 
residual images are compatible with the noise levels of the 
image. These images reveal no evidence for residual emission 
extended beyond the PSF in either band. 

HD 113766 is a known binary star with F3/F5 compo- 
nents. In the Hipparcos catalogue the B component is listed 
at an offset of l'.'335 at a position angle of 281° East of North 
from the primary. A source is clearly seen in the residual im- 
age taken with filter SiC at an offset of l'.'37 at a PA of 279° 
East of North from the primary, based on Gaussian fits to 
determine the centers of the source images. At a parallax 
distance of 131pc this translates to an on-sky separation of 
157AU. Aperture photometry at this location gives a flux 
of 49±5mJy for the binary in the N band. We see no sig- 
nificant flux at this location in the Q band image. This is 
consistent with a SED fit to the binary companion (a Ku- 
rucz profile of spectral type F5 scaled to the Hipparcos V 
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Figure 2. The azimuthally averaged surface brightness profiles 
of the observations of HD111915 (standard) and HD113766. The 
profiles show that the HD113766 images are very similar in shape 
to the standard star, again showing we have no evidence for ex- 
tended emission beyond the PSF. The binary companion surface 
brightness profile is not shown at Q as there is no significant 
emission detected at its location in this band. 



band flux of the binary companion) which predicts a flux of 
14mJy. Such emission is below the la detection threshold 
on the Q band image (Icr in a l"radius aperture is 15mJy). 
The contours on the HD113766B image in the N band (Fig- 
ure[T]top right) are more uneven than the PSF as measured 
on the standard star, but there is no evidence of extended 
emission from the contours. As a further test the PSF image 
scaled to the peak of the binary companion was subtracted 
from the HD113766B image and the residual emission was 
found to be consistent with the noise on the image. 

As a final test for extended emission around HD 113766 
the surface brightness profiles of the standard star images 
and HD 113766 were examined. There is no evidence in these 
profiles for extended emission around HD 113766 beyond the 
shape seen for the point-like standard star targets (Figure 
[2}. If we take sub- integrations of the science observations 
and the standard star observations to consider the variation 
in the profile we have 2 sub-integrations of the standard 
star observation and 5 for HD113766 at N. The mean and 
standard deviations for the FWHM measurements for these 
sub-integrations are 0'.'324±0'.'002 for the standard star and 
0'.'322±0'.'004 for HD113766. For the Q band images we have 
4 sub- integrations of the standard star observation and 11 
for HD113766. The mean and standard deviations of the 
FWHM measurements for these sub-integrations are 0'.'498± 
O'.'OOe for the standard star and 0'.'498±0'.'002 for HD113766. 
These results also show no evidence for extended emission 



Figure 3. The limits on the disc location around HD113766 
placed by non-detection of extension in our images. Different 
source geometries are indicated by different line styles as given in 
the legend. Intermediate inclinations would have limits between 
the face-on and edge-on discs. The regions above the lines rep- 
resent discs that would have been detected at the 3(t level or 
higher. Thus the regions below the lines represent the possible 
disc parameter space, given our non-detection of extended emis- 
sion. The asterisks mark the mid-point o f the p redicted disc lo- 
cation according to the fit bv lLisse et al.l ll2008h . At these levels 
of emission (1599 mjy at N and 1867 mjy at Q) we would have 
expected to detect any disc larger than 0'.'13. 



around HD113766 in either band. The increase in FWHM in 
the Q band compared to the N band is exactly what would 
be expected from the increased diffraction limit at the longer 
wavelength of observation. 

In addition to examination of the images, surface bright- 
ness profiles and FWHM meas u remen ts, we use the tech- 
nique presented in ISmith et al.l (|2008l l to assess the limits 
of our detection capability for extended emission. For each 
photometric band we created a series of disc models com- 
prising rings of radius r with "widths" dr (so that the inner 
radius of the disc would be r — dr/2, outer radius r + dr/2) 
and different inclinations to the line of sight. The rings were 
assumed to have constant surface brightness. These model 
images were added to point sources, representing the stars, 
and these sums then convolved with a PSF model. The fi- 
nal model images of the d iscs plus star s were then subjected 
to the testing detailed in ISmith et al.l l|2008l ). In brief, this 
testing consisted of subtracting the point-like emission (by 
scaling the PSF model to the peak of the image in the same 
way as the residual images above were created), multiplying 
the residual image by a mask, and testing this final image 
for emission above the noise level of the image. The masks 
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Figure 4. The final co-added images for HD171759 (the standard star) and HD172555 observed with TReCS. The contours are at 10, 
25, 50 and 75% of the peak and show that there is no obvious evidence for extended emission beyond the size of the PSF (colour scale 
is linear with brightness). The small area at ~10% of the peak offset from the source in the Q band standard star image is the peak of 
the Airy ring, seen also in the Q band image of HD172555. In the residual image, created by subtracting the final standard star image 
from the HD172555 image after scaling to the peak, there is no evidence for significant residual emission which might be a resolved disc 
in the N band, however we do see lobed emission in the Q band which could be evidence of an edge-on disc. In these images black pixels 
are <la per pixel. 



blank all pixels apart from a region of optimal size and shape 
to detect the disc emission based on the geometry of the disc 
(optimal regions were determined by extensive modelling of 
optimising d isc detectability dep ending on source geometry 
presented in lSmith et al.|[2009bl ). This allowed us to deter- 
mine what disc parameters (and levels of disc flux, which 
were also varied) would have led to a detection of extended 
emission. The resulting 3cr extension limits on the location 
of the excess emission are shown in Figure [31 For a given ra- 
dius the limiting flux required for a 3a detection of extended 
emission is shown for a disc centered on the star with median 
radius at that radius, with a full radial extent dependent on 
the geometry of the disc. The limits suggest that for narrow 
ring models in the N band, with a flux equal to the disc 
flux inferred from the IRS spectrum (1599mJy), any face-on 
or edge-on discs greater than approximately (/.'IS (17AU) 
in radius would have been detected in our data. Similarly 
in the Q band, for discs with 1867mJy of flux (determined 
by subtracting the predicted photospheric emission from the 
Spitzer IRS spectrum of HD113766) any extended disc larger 
than ^tf.'lSS (18AU) in radius would have been visible at the 
3(T level. The photospheric contribution was calculated from 
a scaled Kurucz mode l photosphere as o utlined in Table [T] 
Detailed modelling bv lLisse et al. I (|2008t ) suggests that the 
warm emission comes from a region ~1.8AU (13mas)from 
the primary, and icy grains are situated in a belt at 4-9AU 
(31-69mas) fr om the star . Additional ice at 30-80AU (23- 
61mas) in the lLisse et al.l model contributes to longer wave- 



The science images were also tested in the same way, subtract- 
ing both HD113766A and HD113766B by subtracting scaled stan- 
dard star images. The emission in the masked regions (using a full 
range of masks for a range of disc geometries) was always found 
to be consistent with the noise levels in the images. 



Table 3. Observations imaging HD 172555 and a standard star 
mP 171759) with TReCS, first published in iMoerchen et al.l 

bold) . 



Filter 


Int. time (s) 


Object type 


Star 


Si5 


60 


Cal 


HD171759 


Qa 


60 


Cal 


HD171759 


Qa 


640 


Sci 


HD172555 


Si5 


640 


Sci 


HD172555 


Si5 


60 


Cal 


HD171759 


Qa 


60 


Cal 


HD171759 



length excess but has very low emission at the VISIR wave- 
lengths. The limits from the VISIR imaging p resented here 
are consistent with the fit of iLisse et all 12003). 

The imaging data of HD 113766 would have detected 
extended disc structure at the levels detected in the IRS 
spectrum on scales of >17AU. The lack of evidence for any 
extended emission in either the N or Q band images therefore 
allows us to place an upper limit of 17AU on the extent of 
dust emission around the source. 



2.2 TReCS imaging of HD172555 

We do not have VISIR imaging of HD172555, but re- 
ce ntly 8m-imaging data of this target has been presented 
by iMoerchen et alj (|2010l ). In 680s of on-source integration 
with TReCS in both the Si-5 (Ac = 11.66^m, hereafter N 
band) and Qa (Ac = 18.3^m, hereafter Q band) filters the 
authors found no e viden ce for extended emission. However, 
IPantin fc di Folcd ([2011') recently presented Lucky Imaging 
of this target with VISIR which suggested that HD172555 
appears extended in the Q band, although not in the N 
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Figure 5. Profiles of line cuts 
through individual chop-nod 
integrations in the observa- 
tions of HD172555 and stan- 
dard star HD171759. The left- 
hand column shows line cuts at 
110°EoN, and the right-hand 
column shows cuts perpendicu- 
lar (at 20°EoN) to this. There 
is no difference in the profiles of 
the science and standard star 
target in the N band. In the 
Q band the profiles of the sci- 
ence target show greater vari- 
ation, but overall the profile of 
HD172555 is wider at 110°EoN 
than the standard star tar- 
get. The average profiles shown 
by solid (standard star) and 
dashed lines (HD172555) are 
the mean taken over all indi- 
vidual frames. The FWHM of 
the average profile (determined 
through a 1-d Gaussian fit to 
the profile) is shown by a bar 
in the same linestyle. 



band. To explore this issue f urther, we have obtaine d the 
TReCS raw data presented in iMoerchen et al.l (2010') from 
the Gemini Science Archive. The data were reduced using 
the same custom procedures used for our VISIR imaging of 
HD113766. The science observations were calibrated using 
observations of the standard star HD 171759, taken imme- 
diately before and after the science observations. A sum- 
mary of the observations is given in Table [S] In contrast 
to the VISIR imaging, the TReCS observations were per- 
formed with a parallel chop-nod pattern. As the off-beams 
are unguided which can effect the shape of the PSF, we only 
use the guided images of the targets in our analysis. Sub- 
integrations on the targets were co-added using the Gaussian 
centering techinque employed for the VISIR observations of 
HD113766. Aperture photometry centered on the Gaussian 
peaks of the final images was performed using l"radius aper- 
tures (noise levels were determined from an annulus with in- 
ner radius 2" and outer radius 4"). The calibrated flux found 
in the final images was 1120±67mJy for IID172555 at N and 
1039±85mJy at Q (errors include calibration errors of 3% 
and 8% respectively determined from variation in calibra- 
tion factors between the two standard star observations in 
each filter). The se values are consistent with those quoted in 
iMoerchen et al.l (2010; values were listed as 1155±116mJy 
at N and 1094±164mJy at Q including fiducial 10% and 15% 
calibration uncertainties). 

The final images for the PSF reference (standard star) 
and HD172555 are shown in Figure [H The final images are 
shown with contours at 10, 25, 50 and 75% of the peak. 
There is no evidence for extension in the N band image of 
HD172555. In the Q band image we see some evidence of 
greater ellipticity in the image of HD 172555 than is seen in 
the PSF reference image. This is confirmed in the residual 
image, which is created by subtracting from the science im- 
age the PSF reference image scaled to the peak of the science 
image. The residual image is shown in the right-hand column 



of Figure |4] We see two clear lobes of extended emission in 
the Q band residual image aligned along a position angle of 
110°EoN, but no significant emission in the N band residual 
image. The total fiux subtracted from the Q band image of 
HD 172555 to obtain the residual image is 934 mjy, much 
higher than the predicted flux from the stellar photosphere 
in this fllter of 202mJy. 

To test whether the emission we observed is truly signifl- 
cant, or the result of a varying PSF, we examined the science 
and standard star image profiles in sub-integrations (single 
chop-nod cycles) on the targets. We have 6 sub-integrations 
on the standard star object and 32 on HD172555 in each 
band. Taking a strip centered over the peak of the image of 
3 pixels' width (O'.'09/pixel) and averaging across this width 
we create a 1-d proflle of the image for each sub-integration. 
These proflles are taken at 110° EoN (to coincide with the 
residual emission peaks at Q) and 20° EoN (perpendicular 
to the residual peaks). The proflles are shown in Figure [5l 
We see that in the N band there is no evidence that the pro- 
flles of HD172555 are more extended than the proflles of the 
standard star target. In the Q band the image of HD172555 
is extended in the 110°EoN direction, but not at 20° EoN. 
This suggests that if the extended emission we are view- 
ing arises from a disc, then we are viewing the disc close 
to edge-on, or it is a highly elliptical disc. The proflles in 
the Q band of the science target are quite noisy. As a fl- 
nal test of the extension detection, we examine the FWHMs 
of the sub-integration proflles. The FWHMs are determined 
by fltting a 1-dimensional Gaussian to each proflle. The re- 
sults are shown in Figure [6l It is clear from this plot that 
the FWHMs of the HD172555 images are always larger than 
the FWHM of the standard star images at 110° EoN in the 
Q band. At 20° EoN and at both angles in the N band, the 
FWHMs of HD172555 are consistent with those measured 
for the standard star target. 

The disc models used in the extension limits testing for 
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Figure 6. Measurements of 
the FWHM of the image pro- 
files of HD172555 and standard 
star I^D171759 taken from hne 
cuts through single chop-nod 
integrations centered on the 
stars. The FWHM values were 
determined from a Gaussian 
1-d fit to each profile. Mean 
and rms errors on the FWHM 
values are listed in the leg- 
ends. For the Q band profiles 
at 110°EoN the science target 
HD172555 has a larger FWHM 
in all frames than the stan- 
dard star target. In the N band 
and the Q band at 20°EoN the 
FWHM of HD172555 are con- 
sistent with the FWHMs mea- 
sured for the standard star ob- 
servations. 



HD 113766 were then used to determine the approximate 
size of the extended emission. The models (of varying size, 
thickness, inclination and rotated to different position an- 
gles, see previous section 2.1 for a description) were added 
to point sources representing the star and any unresolved 
excess and convolved with the PSF (the standard star im- 
age). These model images were compared to the image of 
HD172555 by subtracting the model image from the science 
image (scaling to the peak) and comparing the residuals with 
the noise on the image. Using a calculation to determine 
the best fitting model we find that a disc of radius 0'.'27 
(7.9AU at 29.2pc), width dr = 1.2r, inclined at 75° to the 
line of sight lying at a position angle of 120° EoN provides 
the best fit to the data. In this best-fitting model 65% of 
the total fiux arises from the extended emission, suggesting 
a point-like flux of ~ 363mJy, exceeding the 202mJy pre- 
dicted to arise from the star. The best fitting model image 
is shown, after subtracting the scaled PSF image, in Figure 
[7] which can be directly compared to the residual image in 
Figure ID The contours plotted on this figure are from the 
HD172555 residual image, and show that the model does 
indeed reproduce the main features of the extended emis- 
sion. To determine limits on the radius of the disc model, 
we sum the values over all values of the other parameters 
tested {dr, inclination, position angle and flux in the disc). 
Using the percentage points of the distribution we flnd 
that discs with radii Cf.'OQ < r <0'.'31 (2.6-9.1AU) fit the 
observed emission within 3cr. Using the same technique for 
the other model parameters we obtain the following limits: 
0.36r < dr < 2r; inclination > 47°; 40° EoN < position 
angle < 130° EoN; 14% < percentage of total flux in disc < 
74%. 

The focus of this paper is the size and geometry of the 
dust emitting at lO/xm. As we see no extension in the N- 
band image, we use the same procedures adopted in the 
VISIR imaging of HD113766 to place a limit on the size and 
geometry of the HD172555 disc in the N band. The results 
are shown in Figure[8l These limits suggest that the emission 



Figure 7. The best fitting disc model for the Q band extended 
emission seen around HD 172555, shown after subtraction of the 
PSF model for direct comparison to the residual image in Fig- 
ure |4] (colour scale and contour levels the same) . The contours 
overplotted are those from the HD 172555 residual image, and 
show that the residual peaks of the model are of a similar size 
and shape. 



in the N band lies at <0'.'27 (7.9AU) and thus if the mate- 
rial dominating the N band emission was coincident with 
that dominating the Q band extended emission we would 
have expected to detect extension, making this result a sig- 
niflcant non-detection. The size limit matches the size of 
the best fitting model to the extended emission detected in 
the Q band image. As this model showed some evidence for 
residual unextended emission, this may suggest that there 
are two populations of dust in this system. In this case the 
excess around HD 172555 is a possible multiple component 
disc l ike that seen around other main sequence stars (e.g. 
77 Tel lSmith et al.l [2009a ). However, it is also possible that 
we are observing a more extended disc, and that the N and 
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Figure 8. The limits on the disc location around HD172555 in the 
N band based on our analysis of archival TReCS data. Different 
source geometries are indicated by different line styles as given in 
the legend. The regions above the lines represent discs that would 
have been detected at the 3cr level or higher. Thus the regions 
below the lines represent the possible disc parameter space, given 
our non-detection of extended emission. The asterisk marks th e 
predicted disc location according to the fit bv lLisse et al. I l l2009t) . 
At this level of emission (635 mjy at N) we would have expected 
to detect any disc larger than 0''27. The radial offset of the best 
fitting model to the extended emission observed in the Q band is 
shown by a circle. If the material dominating the N band emission 
was coincident with the material which dominates the emission 
at Q we would have expected to detect extension in the N band 
image. 



Q band data are probing the inner and outer parts of the 
distribution due to their greater sensitivity to hotter and 
cooler dust in the disc respectively. This would also account 
for the apparent unresolved component in the Q band im- 
age, as there would still be some emission arising from inner 
part of the disc observable in the Q band. 



2.3 VISIR spectroscopy 

Detailed spectra of both our science targets ha ve been ob- 
taine d with IRS on the Spitzer Space Telesco pe (|Chen et al.l 
bOOeD . As discussed in ISmith et all (2009b) the measure- 
ments of photometry with MIDI suffer from poor back- 
ground subtraction. It is helpful therefore to use the IRS 
spectra as reference total photometry for comparison to the 
correlated fluxes measured with MIDI. However, these spec- 
tra were obtained in a much broader slit than is used in 
MIDI observations (3'.' 7 or 4'.' 7 used in IRS; MIDI slit is 0'.'52 
wide). It is thus possible that extended emission caught in 
the IRS slit would not be observed in the MIDI slit. The 
binary companion to HD 113766 must also be ruled out as a 
source of excess emission. 

To avoid any biasing of the interferometric visibilities 
which could arise from assuming a higher total flux than 
falls within the MIDI slit, we obtained VISIR spectroscopy 
of both science targets in low resolution mode (R~350 at 
lOjum) with a 0'.'75 slit (program ID 083.C-0775(E)). Two 
fllters were used to examine the short and longer wavelength 
ranges covered by MIDI (Ac = 8.8^m range 8-9.6/im, here- 
after fllter 8.8; Ac = 11.4/im range 10. 43-12. 46/im, hereafter 
fllter 11.4). The observations were performed in chop-nod 
mode with standard star observations taken immediately be- 



Table 4. VISIR low-resolution spectroscopy 



Date 


Target 


Filter 


Int. time (s) 


11th May 2009 


HD111915 


8.8 


240 


11^-1^ T\/T„,. onnn 
iitn iViay zOOy 


xlUiiv) (Ob 


8.8 


900 


11th May 2009 


HD111915 


8.8 


240 


11th May 2009 


HD111915 


11.4 


240 


11th May 2009 


HD113766 


11.4 


600 


11th May 2009 


HD111915 


11.4 


240 


25th May 2009 


HD156277 


8.8 


240 


25th May 2009 


HD172555 


8.8 


900 


25th May 2009 


HD156277 


8.8 


240 


7th July 2009 


HD156277 


11.4 


240 


7th July 2009 


HD172555 


11.4 


600 


7th July 2009 


HD156277 


11.4 


240 



fore and after each science observation (Table |4] summarises 
the observations). 

The observations were reduced with 
the VISIR pipeUne procedures available at 
http:/ /www. eso.org/sci/data-processing/software/pipelin es/ 1 
Calibration was performed using an average of the two 
standard star observations obtained either side of the 
scienc e observation s. The standard stars fluxes were taken 
from ICohen et al.l (Il999l ) models for the mid-infrared 
standards used. We plot the observed spectra in Figure [5] 

The VISIR spectroscopy of HD 113766 was performed 
on one night, with the observation block (standard star 
HD111915, science target HD113766 and the standard star) 
in filter 11.4 taken immediately after the observation block 
in filter 8.8. The final observed spectrum looks very simi- 
lar in shape to the IRS spectrum of HD113766, although 
it is lower everywhere by a factor of 1.3. This difference is 
not due to the binary component which fell within the slit 
of the IRS spectroscopy but outside the slit in the VISIR 
observations, because this source contributes only an aver- 
age of 2% to the total fiux across the 8-13/.im range (from 
scaled Kurucz model photos phere as described in Table [T] 
and the previous subsection) . iGeers et al.l l|2007l ) found that 
when using Spitzer IRAC photometry to calibrate VISIR N 
band spectra errors in absolute calibration could be as large 
as 30%, consistent with the difference between the VISIR 
and IRS spectra of HD 113766 seen here. We also obtained 
archival IRS spectra of the standard star HD 111915. As the 
spectral ranges of VISIR and IRS are different, this allowed 
a comparison of the VISIR and IRS spectra of the standard 
star in the longer wavelength observation with VISIR only. 
A comparison of the IRS spectra reduced using pipeline rou- 
tines shows that the calibration factor was varying over the 
course of the longer wavelength observation (we cannot test 
the stability of the calibration for the shorter wavelength ob- 
servations as there is no overlap with the IRS spectra) . The 
difference between the IRS and VISIR spectra was a factor 
of 1.07 and 1.15 for the VISIR calibration observations taken 
before and after the science observation, although again the 
shape of the spectra was consistent between the two. Al- 
though these differences are not as large as the factor of 1.3 
observed for the science target, the difference between the 
two indicates that the absolute calibration was quite unsta- 
ble during these observations. Within the errors of absolute 
calibration for VISIR, we see no evidence for extended emis- 
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Figure 9. The VISIR spectroscopy of HD113766 (left) and HD172555 (right). The IRS spectra of these targets are shown for comparison. 
The stellar emission as calculated from a Kurucz model photosphere of appropriate spectral type scaled to the 2MASS K band emission 
of the star is shown as _Fi. 



sion detected in the IRS spectroscopy that would fall outside 
the MIDI slit. It is also worth noting that we do not see evi- 
dence for temporal evolution in the emission (within the cal- 
ibration errors). The prospect of temporal evolution of the 
excess emission around HD 69830 (another star with bright 
levels of excess in the terrestrial planet region) has recently 
been ruled out (feeichman ct al. 2011), but this remains a 
possibility for a star undergoing planet forming collisions. 
Photometry from th e IRAS database (obtained in 1983; 
iBeichman et al.lll988l ') gives a flux of 1590±80mJy at 12^im. 
The Spitzer Space telescope IRS flux averaged over the fi- 
nite bandwidth of the ^ISIR N band filte r is 1803±92 mjy 
(data obtained in 2004: IChen et al.ll2006l ). The VISIR pho- 
tometry taken just 2 years prior to the VISIR spectroscopy 
presented in section 2.1 (1673±42mJy for HD113766A at N) 
is also consistent with a constant level of excess emission. 

For IID172555 the difference between between the IRS 
spectrum of the target and that measured with VISIR is 
an average of 3% in filter 11.4, consistent with the varia- 
tions observed between different standard star observations 
(measured to be 5%). The observations in filter 8.8 were 
taken on a different night, and the difference between the 
IRS spectrum and VISIR spectrum is found to be closer to 
13%. Taken separately, the two sections of the VISIR spec- 
trum again appear to be simply scaled versions of the IRS 
spectrum. The apparent sharp slightly offset peak in the 
VISIR spectrum at ~9.5/im is consistent with the fiatter 
peak between 9.2-9.5/xm seen in the IRS spectrum within 
the uncertainty on the spectrum. The difference in calibra- 
tion on different nights is the likely cause of the difference 
in scale factors for the observations in filters 8.8 and 11.4. 
In both cases the differences ar e within the 30% a bsolute 
calibration uncertainty found bv lGeers et al.l ( |2007l ). There 
are no IRS observations of the standard star HD 156277 in 
the Spitzer archive. 

Although the errors in absolute calibration for the 
VISIR spectroscopy are large, the shape of the spectra for 
HDl 13766 and HD172555 agree with the shape of the IRS 
spectra for both targets. We would expect cooler emission 
to be further from the star and thus more likely to be ex- 
cluded from the VISIR spectroscopy. As we do not see a 
preferential loss of flux in the VISIR spectra as compared 
to the IRS spectra at long wavelengths, there is no evidence 
on the basis of the spectroscopy for extended emission that 



would fall outside the MIDI slit being detected in the IRS 
spectra. The absolute offsets between the VISIR and IRS 
spectra are consistent with an expected uncertainty of 30%, 
and there is no evidence for temporal evolution in the flux 
levels from near contemporaneous measurements. The IRS 
spectra of both targets shall therefore be used as a mea- 
sure of the expected total photometry for comparison to the 
correlated fluxes measured with MIDI, however the implica- 
tions if the VISIR spectrum had measured the true flux of 
the HD 113766 system are discussed briefly in sections 3.2 
and 4.3.1. 



3 THE MIDI OBSERVATIONS 

Observations were taken over several semesters through a 
combination of service and visitor mode observations. Table 
[5] lists the observing run IDs and dates of all observations, 
together with the baseline conflgurations used. All obser- 
vations used MIDI on the UTs in the HIGH-SENS mode 
where the interferometric fringe exposures are followed by 
separate photometric exposures from each telescope in turn 
to measure the target spectrum thr ough each beamline (for 
a summary see ISmith et al.l l2009bl : furthe r details can b e 
found in the MIDI instruction manual or in iTristramllioOTi '). 

Reduction of the MIDI data was performed using the 
EWS software, available as part of the MIA+EWS package 
(see http:/ /www. strw.leiden.nl/~nevec/MIDI/index. html). 
Reduction followed the standard EWS routines. A summary 
of these steps is given below. Further details can be found 
in the manual available at the above link. 

• Frames were multiplied by a mask and compressed in 
the direction perpendicular to the spectral dispersion to ob- 
tain a one-d i mensio nal fringe intensity spectrum. Following 
ISmith et all (|2009bl ) we used masks determined from a fit to 
the total intensity (photometry) frames, as provided by the 
MIA reduction package, to best exclude source-free back- 
ground pixels and thereby reduce noise levels. 

• The interferometric fringe data were aligned in time 
using an analysis of the measured group delay to remove 
components of the delay arising from instrumental and at- 
mospheric effects. For further details of this procedure the 
reader should consult the EWS manual or Tristram (,2007l ). 
Fringes were then averaged in time to produce a correlated 
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Table 5. Observations of science and calibrator targets with MIDI 



Date 


Observing 


Baseline 


Baseline 


Baseline 


Target 


Target 


Seeing 


TO 


Flux 




ID 


configurations 


length (m) 


Pos. Angle (°) 


name 


type 




ms 


RMS 


08/03/2007 


078.D-0808(D) 


UT1-UT3 


102.2 


36.52 


HD116870 


Cal 


1.76 


2.0 


0.0085 


08/03/2007 


078.D-0808(D) 


UT1-UT3 


92.8 


0.97 


HD172555 


Sci 


0.70 


6.2 


0.0034 


09/04/2007 


079.C-0259(G) 


UT1-UT3 


101.8 


2.8 


HD 113766 


Sci 


0.72 


3.1 


0.0024 


09/04/2007 


079.C-0259(G) 


UT1-UT3 


102.2 


9.0 


HD112213 


Cal 


0.83 


2.6 


0.0053 


10/04/2007 


079.C-0259{G) 


UT1-UT3 


101.8 


-2.7 


HD113766 


Sci 


0.59 


3.0 


0.0362 


10/04/2007 


079.C-0259(G) 


UT1-UT3 


102.3 


3.0 


HD112213 


Cal 


0.74 


2.8 


0.0028 


30/05/2007 


079.C-0259{F) 


UT1-UT2 


48.7 


41.7 


HD112213 


Cal 


1.52 


1.0 


0.0129 


30/05/2007 


079.C-0259(F) 


UT1-UT2 


42.6 


49.0 


HD113766 


Sci 


1.73 


0.7 


0.0297 


18/03/2008 


080.C-0737(C) 


UT3-UT4 


57.9 


85.8 


HD156277 


Cal 


0.81 


5.7 


0.0034 


18/03/2008 


080.C-0737(C) 


UT3-UT4 


56.1 


78.8 


HD172555 


Sci 


0.59 


7.6 


0.0020 


20/03/2008 


080.C-0373(D) 


UT1-UT3 


88.9 


20.4 


HD156277 


Cal 


0.88 


6.4 


0.0029 


20/03/2008 


080.C-0373(D) 


UT1-UT3 


92.6 


7.8 


HD172555 


Sci 


0.61 


9.7 


0.0021 


20/03/2008 


080.C-0373(D) 


UT1-UT3 


100.7 


12.8 


HD169767 


Cal 


0.51 


9.0 


0.0026 


20/03/2008 


080.C-0373(D) 


UT1-UT3 


92.0 


14.5 


HD172555 


Sci 


0.49 


10.3 


0.0025 


21/03/2008 


080.C-0737(E) 


UT2-UT4 


89.3 


60.8 


HD156277 


Cal 


0.82 


4.6 


0.0029 


21/03/2008 


080.C-0737(E) 


UT2-UT4 


88.1 


45.6 


HD172555 


Sci 


0.71 


4.8 


0.0034 


21/03/2008 


080.C-0737(E) 


UT2-UT4 


89.4 


70.7 


HD156277 


Cal 


0.79 


5.8 


0.0023 


07/05/2009 


083.C-0775(C) 


UT3-UT4 


58.7 


90.0 


HD171759 


Cal 


1.29 


1.8 


0.0037 


07/05/2009 


083.C-0775{C) 


UT3-UT4 


59.3 


95.5 


HD172555 


Sci 


1.10 


1.8 


0.0044 


07/05/2009 


083.C-0775(C) 


UT3-UT4 


60.8 


103.9 


HD171212 


Cal 


1.10 


1.7 


0.0040 


07/05/2009 


083.C-0775(C) 


UT3-UT4 


60.2 


102.7 


HD171759 


Cal 


1.14 


2.0 


0.0036 


07/05/2009 


083.C-0775{C) 


UT3-UT4 


60.9 


108.0 


HD172555 


Sci 


0.95 


2.0 


0.0044 


07/05/2009 


083.C-0775{C) 


UT3-UT4 


62.4 


135.4 


HD152186 


Cal 


0.82 


2.0 


0.0035 


07/05/2009 


083.C-0775(C) 


UT3-UT4 


61.2 


114.7 


HD171759 


Cal 


0.98 


2.3 


0.0027 


07/05/2009 


083.C-0775(C) 


UT3-UT4 


62.1 


121.4 


HD172555 


Sci 


1.49 


1.5 


0.0039 


07/05/2009 


083.C-0775(C) 


UT3-UT4 


62.4 


128.2 


HD171212 


Cal 


1.14 


1.9 


0.0024 


07/05/2009 


083.C-0775(C) 


UT3-UT4 


62.0 


130.2 


HD171759 


Cal 


0.76 


2.9 


0.0026 


08/05/2009 


083.C-0775(D) 


UT1-UT3 


86.2 


49.6 


HD112213 


Cal 


0.47 


6.2 


0.0017 


08/05/2009 


083.C-0775(D) 


UT1-UT3 


82.2 


52.6 


HD113766 


Sci 


0.46 


6.7 


0.0026 


08/05/2009 


083.C-0775(D) 


UT1-UT3 


75.8 


55.7 


HD110253 


Cal 


0.52 


6.8 


0.0015 


08/05/2009 


083.C-0775(D) 


UT1-UT3 


75.7 


55.2 


HD112213 


Cal 


0.51 


5.8 


0.0021 


08/05/2009 


083.C-0775{D) 


UT1-UT3 


72.9 


57.9 


HD113766 


Sci 


0.49 


6.2 


0.0020 


09/05/2009 


083.C-0775(B) 


UT2-UT4 


74.6 


113.0 


HD112213 


Cal 


0.97 


3.0 


0.0027 


09/05/2009 


083.C-0775{B) 


UT2-UT4 


74.9 


116.2 


HD 113766 


Sci 


1.16 


2.6 


0.0020 


09/05/2009 


083.C-0775(B) 


UT2-UT4 


67.6 


126.0 


HD110253 


Cal 


0.71 


4.9 


0.0022 


09/05/2009 


083.C-0775(B) 


UT2-UT4 


69.9 


125.2 


HD 113766 


Sci 


0.76 


4.9 


0.0022 


09/05/2009 


083.C-0775{B) 


UT2-UT4 


62.6 


132.9 


HD112213 


Cal 


0.71 


4.8 


0.0020 


09/05/2009 


083.C-0775{B) 


UT2-UT4 


89.4 


65.9 


HD171759 


Cal 


0.56 


5.9 


0.0020 


09/05/2009 


083.C-0775(B) 


UT2-UT4 


89.4 


68.8 


HD172555 


Sci 


0.62 


5.6 


0.0025 


09/05/2009 


083.C-0775(B) 


UT2-UT4 


89.4 


74.6 


HD171212 


Cal 


0.80 


4.5 


0.0030 


09/05/2009 


083.C-0775(B) 


UT2-UT4 


89.2 


78.0 


HD171759 


Cal 


0.56 


5.6 


0.0034 



Seeing, tq (coherence time) and flux RMS are taken from the ESO ambient conditions database 

(http:/ /archive. eso.org/cms/eso-data/ambient-conditions ) for the time at which the interferometric stage of the observations was taken. 
Seeing is defined as the FWHM of a stellar image observed with an infinitely large telescope at 500nm wavelength and at the zenith. 
Flux RMS gives a measure of the background, with levels >0.05 indicating cloud cover and levels >0.02 indicating possible cloud cover. 
Coherence times (ro) of less than 3ms are considered very fast, and indicate the presence of rapid atmospheric fluctuations that are 
likely to have degraded the interferometric signal-to-noise. 

Science observations were calibrated using the bright standard star observations obtained before and after the science observation 
where possible (i.e. when the science observation occured between two standard star observations close both in time and on-sky 
position). Otherwise only the observation of a bright standard star (standard taken from the CalVin tool, see Table [T]| closest in time 
to the science observation was used for calibration purposes. 
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Figure 10. The visibility transfer function as measured on bright standard star targets cahbrated by other bright standard star 
observations. The region between 9.2-10/im is subject to high ozone absorption and is therefore excluded from the analysis. The dashed 
lines represents an observation of HD112213 taken under proposal 079.C-0259{G) calibrated by an observation of the same target taken on 
the same baseline configuration but on the following night. The dotted line shows an observation of HD112213 observed on 083.C-0775(B) 
calibrated by an observation of HD171759, which is > 90° away, but taken immediately following on the same baseline configuration. In 
the right hand panel the transfer functions are scaled to 1 in the wavelength range 8-9.2/im to show how the transfer function changes 
with wavelength. These observations are discussed in the text. 
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Figure 11. The visibility transfer function as measured on faint standard star targets calibrated by bright standard star observations. 
The presentation of the data is the same as the previous figure. For these fainter targets, the transfer function is flat between 8.5 and 12 
microns, but then shows a systematic increase at longer wavelengths, where in addition the signal-to-noise is starting to degrade. 



flux (or more correctly correlated intensity Zcorr as no flux 
calibration had been determined at this point). The corre- 
lated flux Fcoii was then compared to the total source flux 
Ftot to give the source visibility V = -FcorrZ-Ftot- 

In principle, the total source flux (Ftot) could have been 
determined from the photometry frames observed following 
the fringe exposures with MIDI. The MIDI photometric data 
are consistent with the IRS spectra, however the variation 
between individual observations of the science target pho- 
tometry with MIDI was as high as 30-40% across the full 
wavelength range. We have found that photometric mea- 
surements of faint targets wi th MIDI can often be rather 
noisy (see ISmith et aLlbOOgbl ) and so instead we used the 
IRS spectra — which are not significantly different from the 
VISIR measurements (see Section 2) — to provide these 
data. 

3.1 Visibility calibration 

Two types of standard stars were used to calibrate 
the observations of the science targets. Bright stan- 
dards were selected with the ESO CalVin tool (see 



"http:// www.eso.org/instruments /midi / tools ) . Much fainter 
standards were identified by searching for sources in the 
IRAS catalogue, within 25° of the science targets, that had 
similar 12^m fiuxes and that showed no evidence of bina- 
rity or variability. Our decision to utilise both faint and 
bright calibrators was motivated by a finding from our ear - 
lier studies of HD69830 and rj Corvi l|Smith et al.l l2009bl ) 
which showed tentative evidence of a loss of correlated flux 
at shorter wavelengths for faint targets. Our goal was to 
observe faint standards as though they were famt science 
targets to check for any bias in their measured visibilities. 
For these tests, we ensured that the 12/im fluxes of the faint 
stand ards were con s istent with the emission predicted from 
scaled ICohen et al.1 l|l999t ) photospheres, and were confident 
that they showed no evidence of any excess emission (see 
next paragraph). 

For those targets in the lCohen eT^ (|l999D catalogue of 
mid-infrared standard stars (HDs 111915, 116870, 156277, 
112213) we used the Cohen spectrum of the target as the tar- 
get fiux. For the remaining targets we used Cohen templates 
for stars of the same spectral type scaled to the lO/^m fiux 
listed in the CalVin tool (for HDs 169767 and 171759) or for 
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the faint standard stars (HDs 171212, 152186 and 110253) 
scaled to the source's listed 2MASS K band flux. The di- 
ameters for the bright standard stars used for the MIDI ob- 
servations were taken directly from the CalVin tool. For the 
faint standards source diameters were determined by assum- 
ing that the stars had a diameter typical for their spectral 
type (taken from ICo^ lioOO) and using the Hipparcos listed 
parallax to determine their distances. Our inferred source 
diameters for the standards are listed in Table [T] 

The calibration of the target correlated flux was deter- 
mined using the following equation: 

-Pcorr,tar ~ (-^corr,tar /-^corr ,cal ) ^ -^tot,cal ^ ^al 5 (1) 

where Icon, tar and /corr,cai werc the correlated intensities 
measured in the fringe exposures of the 'target' and 'calibra- 
tor' respectively and Vcai was the visibility of the calibrator, 
assumed to be that of a uniform disc with the diameter given 
in Table [T] For Ftot,cai we used the total flux of the calibra- 
tor taken from the Cohen spectrum of the standard (or the 
scaled Cohen spectrum, see paragraph above) . The visibility 
of the target was then evaluated as 

Vtar — -fcorr,tar/-ftot,tar 5 (2) 

where -Ftot.tar was the total flux of the target, which in the 
case of the science targets was taken from the Spitzer IRS 
spectrum of the source. 

To determine the accuracy of our derived visibilities, 
we first examined the visibilities of the bright standard star 
targets when calibrated by other bright standard stars. By 
using pairs of standard star observations taken as close as 
possible in time, we were able to generate 10 independent 
visibility functions. These pairs of standard star observations 
were in general taken either side of science observations, and 
so were separated in time by roughly twice the time between 
a science and standard star observation. When corrected for 
the sizes of the standard stars used, these produced the vis- 
ibility transfer functions shown in Figure [TOl Note that the 
bandpass between 9.2-10/im is subject to high levels of un- 
certainty due to ozone absorption, and so this region has 
been ignored in our analysis. The variation in transfer func- 
tion level in the left hand panel of Figure [10] is indicative 
of the range of seeing mismatch between the observations of 
source and calibrator, with the values nearest to unity being 
associated with the use of calibrator stars closest in time and 
space. The weighted mean of the transfer functions across 
the whole spectral range gave a mean value of observed visi- 
bility/predicted visibility of 1.002±0.057 (or ~ 6%). To cal- 
culate the weighted mean we used weights derived from the 
errors on the correlated flux. These errors came from the 
variance found in the EWS reduction by splitting the fringe 
observations into 5 sub-integrations. We then calculated the 
weighted mean over the MIDI spectral range (excluding the 
9.2-10/im region which suffers from ozone absorption), and 
took a mean over all 10 pairs of standard stars to get the 
figure above. 

In order to remove this source of variation, we scaled 
each transfer function so that its weighted mean in the range 
8-9/im was unity. These normalised transfer functions are 
presented in the right-hand panel of Figure 1101 Most of the 
normalised transfer functions are very similar, but two of 
them deserve mention. The first, identified by a dashed line, 
is from an observation of HD112213 taken on 09/04/2007 



(under observing ID 079.C-0259(G)) calibrated by an obser- 
vation of the same target taken under the same observing 
configuration the next night. The second was derived from 
an observation of HD112213 taken on 09/05/2009 (under 
proposal 083.C-0775(B)) but has been calibrated using a 
"standard" located over 90° away on the sky (HD 171759). 
These aberrant transfer functions highlight the need for very 
careful calibration strategies. Overall, we found the ratio of 
the transfer function in the range 10.5-11.5/im (12-13/im) 
to that at 8-9^m for the remaining bright-bright pairings to 
be 1.010±0.013 (1.016±0.018). These data confirm that for 
our data the visibility functions for bright targets are likely 
to be calibrated to within ~6%, and also that the differen- 
tial visibility (the visibility with reference to that at a fixed 
wavelength) is a factor of 2-3 times more accurate than the 
absolute value of the visibility. 

A similar analysis for our faint unresolved targets cali- 
brated by bright standard stars is shown in Figure [TTI The 
behaviour of these transfer functions is broadly similar, but 
there is definite evidence of curvature of the functions at 
the extremes of the MIDI bandpass. This is limited to the 
very first few spectral bins at the short wavelength end, 
but is more noticeable beyond a wavelength of 12/xm, where 
the data are increasingly noisy. We found that the weighted 
mean absolute values of the transfer functions across the 
MIDI wavelength range (excluding the 9.2-10/.tm region) 
was 1.035±0.072, and that the mean normalised transfer 
function between 10.5-11.5/xm relative to that at 8-9/im was 
0.984±0.030. The equivalent value for the U-lSfim band- 
pass was 1.064±0.047. These data suggest that the visibil- 
ities of our faint targets can be measured to better than 
10% (i.e. at a level c onsistent with what has been presented 
before, see e.g. [Ch esnoau 2007) , and that the differential vis- 
ibility across the whole wavelength range is stable to better 
than 5%. Somewhat better accuracy can be expected if data 
at wavelengths >12^m is excluded, typically by a factor of 
two. 

3.2 MIDI observations of HD 113766 

The observations of HD 113766 were calibrated using two 
bright standard star observations where possible (see Table 
[5ll. In these instances the average of the calibrated corre- 
lated flux using two standard star observations was used. 
The standard deviation between the two calibrations was 
added in quadrature to the error on the correlated flux as 
estimated from the time variation in the correlated inten- 
sity. This was derived by splitting the fringe integration into 
5 sub-integrations and determining the standard deviation 
between the different sub-integration datasets. 

The correlated ffux measurements for our 7 observations 
of HD113766 taken on various baselines (see Table [5} are 
shown in Figure [T2l where the solid line is the IRS spectrum 
of the target. It is immediately clear that there is a signifl- 
cant difference between the correlated flux and the IRS spec- 
trum. There also appears to be some difference between the 
observations taken under proposal 079.C-0259 and proposal 
083.C-0775. This is further explored through comparison of 
their visibility functions. 

The visibilities calculated from the correlated flux mea- 
surement are shown in Figure [121 Interestingly, the visibility 
on baseline 79-F is somewhat lower than that measured on 
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Figure 12. The correlated flux measurements (top left) and subsequent visibility curves (top right) for the MIDI observations of 
HD113766. In the left hand top panel the IRS spectrum (thick black line) is compared to the correlated flux measurements for reference. 
For the visibility curves the Vtbb prediction assumes that all the excess emission detected has been resolved, and that the only contribution 
to the correlated flux would be emission from the star itself. This is not a flat line as the relative contribution of star and excess emission 
to the total flux from the system changes through the MIDI wavelength range (see Figure |9ll. The maximum predicted visibility (a flat 
line at 1) assumes none of the emission has been resolved. The labels for the different observations are coded according to the semester 
during which the observations were taken, the run letter and a number indicating which observation on a particular run it is, so for 
example the flrst observation of HD113766 taken under proposal 079.C-0259(G) is given the label 79-Gl. The baselines of all observations 
of HD113766 are shown in the bottom panel. For the observational results the same colours and line-styles are used for all three plots. 



the longer but approximately parallel baseline 83-D (see cap- 
tion to Figure [121 for description of labels and Table [5}. This 
might indicate that rather than having a distribution like a 
Gaussian, the source emission is more complex. A ring for 
example would have an oscillating visibility function which 
would be higher on some longer baselines (see, e.g.. Figure 5 
of lPuUem ond fc Monn icr 201Q for an example of this) . How- 
ever, the observations on 79-F were taken under the poorest 
observing conditions in the study — the "flux RMS" and 
coherence time diagnostic metrics were particularly high — 
and under these conditions it is likely that the measured 
visibility was biased to a lower value. As a result, we did 
not use this visibility measurement as a constraint in the 
modelling of the HD 113766 data presented in Section 4. 

In the plot of the visibility function (Figure[T2]top right) 
we have shown the behaviour for two extremal classes of tar- 
gets, i.e. those that are completely unresolved and then those 
whose dust emission is fully resolved. For the first of these 
cases, we expect that the correlated fiux is always equal 
to the total fiux, and hence that the visibility is equal to 
unity at all wavelengths. In the second case, the situation 
is slightly more complicated, since both the potentially un- 



resolved stellar contribution and the resolved dust emission 
need to be considered. For such targets, the visibility of the 
source will be given by 

V = I^V. + §if^l/di.c, (3) 

where Vdisc and V* are the visibility of the disc and the star 
respectively, F^isc and F^, are the fluxes of the disc and star 
respectively, and Ftot is the total flux (i.e. Ftot = ^disc + J*)- 
The stellar emission component in HD113766 is likely to be 
completely unresolved. It is a F3V-type star at a Hipparcos 
distance of 131pc, and so with an expected radius of 1.38i?0 
would subtend an angle of roughly 0.048mas, well beyond 
the resolving power of VLTI/MIDI. The minimum expected 
visibility, given that the star is expected to be completely 
unresolved (V* — 1), is therefore the visibility in the event 
that the disc flux is completely resolved, Vdisc = 0, V = 
Vies = -F*/Jtot- We have labelled this visibility Kcs in Figure 
1121 As it is clear that for all baselines observed the visibility 
function lies between this and the maximum visibility of 1, 
the disc appears partially resolved on all baselines. 

The visibility functions measured on similar baselines 
are consistent within the error levels expected for the visi- 
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Figure 13. The weighted mean of the visibility functions (over the MIDI wavelength range) measured for HD113766 shown against 
baseline length (left) and position angle (right). Means are weighted by errors added in quadrature from the statistical error on the 
correlated flux, error on the calibrator correlated flux, error on the calibrator total flux and error on the IRS spectrum used to derive the 
visibility from the correlated flux (see Section 4.2 for more details). Symbol colours are given in the key of Figure [121 The observation on 
baseline 79-F is excluded from this figure due to the high uncertainty associated with this data point (see text for details) . The dashed 
line marks the visibility expected for a model in which the excess emission is distributed as a circularly symmetric Gaussian with FWHM 
lOmas. This gives a first order approximation for the size of the emitting region. Also plotted in grey dashed lines are 5o" limits to the 
Gaussian model (see section 4.3). The best fitting ring-like model for the distribution is also shown by a dotted line. See section 4.3 for 
a detailed discussion of these fits to the observed visibilities. 



bility of a faint target (~10%, see Section 3.1), The visibility 
functions measured on 79-Gl and 79-G2 are very similar to 
one another, as are the visibility functions measured on 83- 
Bl and 83-B2. It is also clear that the visibility functions 
measured on longer baselines are lower than those measured 
on shorter baselines (with the exception of 79-F as discussed 
above, and which we treat as unreliable). Solid lines are used 
in Figure [T2I to denote the longest baselines, dotted lines the 
shortest and dashed lines intermediate baseline lengths. 

To display more clearly how the visibilities are varying 
with baseline length and position angle we show the average 
(weighted mean over the MIDI spectral range excluding the 
ozone dominated 9.2-10/^m region) visibility of HD113766 
plotted against baseline length and position angle in Fig- 
ure [TS] Although the visibility is seen to change with wave- 
length (Figure [121 top right), to first order this effect can be 
attributed to the decreasing resolution of MIDI with increas- 
ing wavelength, and so we are unlikely to be biasing the data 
through this averaging. Excluding baseline 79-F, there is a 
drop in visibility with increasing baseline length, consistent 
with a simple source geometry such as a Gaussian. Also 
shown in this figure are the visibilities plotted as a function 
of baseline position angle. Although the visibilities closest to 
or 180° appear lower, these are the longest observed base- 
lines. If we consider only baselines of similar lengths (83-Dl, 
83-D2, 83-Bl and 83-B2) there is no evidence for a change 
in visibility with baseline position angle that would indicate 
a non-circularly symmetric source distribution. More details 
on the limits we can place on the source geometry with these 
observations are discussed in Section 4. 



^ If we had instead used the VISIR spectrum as a model for the 
total flux of the HD 113766 system, then the calculated visibilities 
would be higher, but still exhibit the same behavior, i.e. decreas- 
ing with baseline length. In this case the points would lie close to 
the 7mas model shown in Figure [T3l 



3.3 MIDI observations of HD172555 

Our observations of HD 172555 have been taken over several 
semesters. Where possible, two bright standard star obser- 
vations were used to calibrate the correlated flux measure- 
ments of this target (see Table [5| ■ Errors were calculated 
in the same way as described for HD 113766 in the previous 
subsection. 

The calibrated correlated flux measurements from our 
9 observations of HD 172555 with MIDI (see Table ^ are 
shown in Figure [141 As for HD113766, the IRS spectrum has 
been over-plotted for comparison. These plots show rather 
similar correlated flux levels for all baselines. The degree of 
resolution of the target is better seen in the visibility func- 
tions of the target, also shown in Figure [141 in the right hand 
panel. As for HD113766, we have over-plotted the predicted 
visibilities for targets with fully unresolved {V = 1 across 
all wavelengths) and fully resolved (Kcs) disks for reference. 

These plots show that for observations 78-D and 83-C3 
the observed visibilities are more than 10% lower than the 
predicted visibility in the case that the disc is fully resolved 
(Vrcs). Although values lower than Kes can be observed on 
intermediate baselines (when Vdisc < due to the phase 
jump s in its pattern, see Figure 5 of lDuUemond fc Monnien 
I2OI0I for an example of this behavior), the visibilities ob- 
served on similar baselines (e.g. 81-Dl and 83-C2 respec- 
tively) do not fall significantly below Kcs. We believe that 
the results on baselines 78-D and 83-C3 can be explained 
by a combination of calibration errors. The coherence times 
for the observations on baseline 83-C were very short, par- 
ticularly for the third observation of HD 172555 which has 
a significantly shorter coherence time than the observations 
of standard stars used to calibrate the observation (see Ta- 
ble O . This could be the cause of the low visibility as we 
believe was the case for our observation on baseline 79-F 
for HD113766. Secondly, the observation on baseline 78-D 
was taken at the end of a night towards twilight and had to 
be calibrated by observations of a standard star offset from 
HD172555 by 95°. As demonstrated in Section 3.1 large off- 
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Figure 14. The correlated flux measurements (top left) and subsequent visibility curves (top right) for the MIDI observations of 
HD172555. In the top left panel the IRS spectrum (thick black line) is compared to the correlated flux measurements for reference. For 
the visibility curves the Vi-es prediction assumes that all the excess emission detected has been resolved, and that the only contribution 
to the correlated flux would be emission from the star itself. This is not a flat line as the relative contributions of the star and excess 
emission to the total emission vary with wavelength (see Figure [9||. The maximum predicted visibility (a flat line at 1) assumes that 
none of the emission has been resolved. Different runs are labelled following the convention described in Figure [12] The baselines of all 
observations of HD172555 are shown in the bottom panel. For the observational results the same colours and line-styles are used for all 
three plots. 
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Figure 15. The weighted mean of the visibility functions (over the MIDI wavelength range) measured for HD172555 shown against 
baseline length (left) and position angle (right). Means are weighted by errors added in quadrature from the statistical error on the 
correlated flux, error on the calibrator correlated flux, error on the calibrator total flux and error on the IRS spectrum used to derive 
the visibility from the correlated flux (see Section 4.2 for more details). Symbol colours are given in the key of Figure [Ml The data from 
baselines 78-D and 83-C3 are excluded from this plot due to high uncertainty associated with these observations (see text for details). 
The dashed line indicates the mean visibility (across the MIDI wavelength range) in the case that the excess emission is completely 
resolved (Vros), indicating that a model in which the disc has a large spatial extent is likely to provide a good flt to the data. There is 
tentative evidence (at the 1-2(7 level) that the disc emission is not fully resolved on baselines 80-D2 and 80-E, suggesting the emission 
may not be circularly symmetric. 
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sets between a science target and standard star are an obvi- 
ous source of miscalibration, which we believe has occurred 
here. 

We can see that for all our baselines the observed visi- 
bility function for HD 172555 is consistent with a completely 
resolved disc. Ifigh visibilities at >11.5/im are likely due to 
the bias seen in observations of faint targets and can be 
compared with those seen for the faint calibrators (Figure 
[TT] right). There is no evidence of a systematic decrease in 
visibility with increasing baseline length, suggesting that the 
excess emission is already completely resolved on the short- 
est baselines. There is tentative evidence that the mean ob- 
served visibility changes slightly with baseline position an- 
gle, which can be seen more clearly in Figure 1151 Such a 
change with baseline position angle could reveal evidence of 
a clumpy structure, as may be expected as the result of a 
recent mas si ve co llision in the disc, the model favoured by 
iLisse et all (|2009l ). However, as the visibility seems to be 
smoothly changing with position angle (increasing with in- 
creasing position angle from 0-14° on baselines of similar 
length 80-Dl and 80-D2, and decreasing with increasing po- 
sition angle from 46° on baselines of similar length 80-E and 
83-B) this suggests a smoother, perhaps, elliptical structure 
for the emission. An ellipse with major axis oriented at 120° 
would have its lowest visibility on baselines at 120° EoN and 
its highest visibility on baselines at 30° EoN. Such an emis- 
sion morphology could arise from an inclined circular disc or 
a truly elliptical disc. As we have complete resolution of the 
disc on the shortest baselines (56m, 80-C) we would expect 
such a disc to have a minimum radial size of at least ~40mas 
(based on the resolving power of 56m aperture telescope at 
10.5/im). More detailed modelling of the observed visibility 
functions is presented in the following section. 



4 MODELLING THE VISIBILITY FUNCTIONS 
4.1 Visibility calculation 

We consider several source geometries to try to fit the ob- 
served visibility functions around HD113766 and HD172555. 
The van-Cittert Vernicke theorem states that the normalised 
visibility function of a source is the normalised Fourier trans- 
form of the brightness distribution of the source. The sim- 
plest source geometry we consider is a circularly symmetric 
Gaussian. The source geometry is given by 



/(q,/3) = 



1 



■ exp 



-4 In 2p^ 

—Q2 



(4) 



^7r/(41n2)e 

where ct and /3 are angular on-sky coordinates (in radi- 
ans), O is the full- width at half- maximum (FWHM) and 
p — yQ^ + /3^. The visibility function of this source is then 



V{u, v) — exp 



41n2 



(5) 



where u and v are coordinates describing the spatial fre- 
quency of the brightness distribution such that u = Bu/\, 
V = By/\ where Bu and B^ are the projections of the base- 
line vecto r on the two axes and A is t he wavelength of the ob- 
servation (|Berger fc Segransanll200"7t ). For an elliptical Gaus- 
sian we parameterise the ellipticity by I where vi = vcosi 
(and ui — u). 



The Gaussian model provides a good model to an enve- 
lope, and offers a simple approximation to the overall source 
size. However, debris disc emission is expected to be ring- 
like in distribution. For a thin ring the source geometry is 
given by 



Hp)- 



1 



5{p- po) 



(6) 



(where 5 is the Dirac delta function), then the visibility of 
such an object is given by 



V{u,v) = Jo{2-Rpor) 



(7) 



where Jo is the Oth-order Bessel function and r — yjv? + v'^. 
For our deb r is dis c models we follow the example of 
iMalbet et af] l|2005h and use a sum of thin rings to model 
the distribution and visibility function of a ring of finite 
thickness. The emission is distributed as the integration of 
all ring contributions given by equation |6] from po = Pmin to 
po = Pmax, and the visibility function of this finite thickness 
ring is similarly the integral of the corresponding Fourier 
transforms (visibility functions). For a disc inclined to the 
line of sight at an angle I and at a position angle 9 we sim- 
ply consider the case where rgj = y^Ug + Vg cos'^{I) which 
represents the projected baseline in a new {ue,vg) reference 
frame corresponding to a rotation of the array frame by the 
position angle 9 with a compression factor of cos(7). To in- 
clude the temperature distribution of the disc, we make the 
simplifying assumption that the disc is optically thin and 
the dust behaves like a blackbody, and so the temperature 
of each ring is dependent only on the distance from the star 
(as the stellar luminosity is fixed). Then 



T(r) = 278.3 



(8) 



and the flux we can expect from our disc model is the inte- 
gration of all the rings modified by the distance of the source 
from the observer, d, 



Fx 



2n 
d2" 



rBA(T(r))Eodr. 



(9) 



We assume the disc surface density is flat and therefore Eq 
is a constant. The visibility function corresponding to this 
model is therefore simply the normalised integration of the 
visibility function for the thin rings weighted by their flux. 

These models are simple descriptions for the visibility 
function of the excess emission, Vdisc- In our comparison to 
the data the final visibility model includes the contribution 
from the star which is completely unresolved by the inter- 
ferometer (V* ~ 1, see Section 3.2). The final value of Vmod 
is c alculated ac c ording to equation O 

iLisse et al.1 (|2009l ') suggested that a possible origin for 
the emission around HD172555 is a recent massive collision 
between two large planetesimals/proto-planets. In this sce- 
nario we might expect the emission to arise from a clump 
offset from the central star. As the clump is offset from the 
central star, we cannot simply add the visibility functions 
scaled by their relative fiux levels. For a multi-component 
function components at positions a^, in the plane of the 
sky with visibilities Vi, the normalised visibility of the full 
function is 



V(u,v) = 



(10) 
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where Fi is the flux of the component i and n is the total 
number of components. For a two-component model, the 
normalised squared visibility reduces to 

FM + FiVi +2F^F2ViV2Cos.{2Tx{u[ar~a2)+v[pi-p2))) 
(Fi + F2Y 

which reduces to the familiar (Fi Vi + ^2^2)^/(^1 + ^2)^ in 
the case that cti = 02 and /3i =132- 



4.2 Error calculation 

To determine the goodness of fit of the models, and thereby 
calculate the best fitting model parameters, the calculation 
of the error terms on the visibility must be carefully con- 
sidered. From equation [T] the error on the correlated flux 
arises from the terms /corr.tar, /corr.cai and -Ftot.cai- Here we 
have assumed that the error on the visibility of the calibra- 
tor is neglig ible. Errors on the calibrator diameters can be 
up to ±5% ( Vcrhoclst"2005') , and so at its highest the error 
on Veal (from observation of HD112213 on baseline 79-G) 
is 0.5%, much lower than the other sources of uncertainty. 
The uncertainty on the cali brator flux is assum ed to be 2% 
(5Ftot,cai/Kot,cai = 0.02, scc ICoheu et"ZI |l999i). 

The error on 7corr,tar is determined by splitting the 
fringe observation of the target into 5 sub-integrations and 
determining the standard deviation from the mean value at 
each wavelength. The error on Icorr.cai is the error from cali- 
bration. Here we use the visibility transfer functions (Section 
3.1) to determine the error. At each wavelength value sam- 
pled in the MIDI range, we determine the mean and stan- 
dard deviation of the visibility transfer functions for the faint 
standard stars. We check that the mean value is compati- 
ble with a value of 1 within the error given by the standard 
deviation, and adopt this standard deviation as our error 
on /corr,cai- The average value of this error across the MIDI 
range (excluding the ozone dominated 9.2-10/xm region) is 
0.057. The error on the correlated flux is then given by 



paragraph. The errors on F-t and Fdisc are more complex. 
They arise in part from the uncertainty on the total flux, 
and also from uncertainty in the relative contributions of 
the star and disc to the total flux. As a change in the value 
(11) of Ftot,tar would result in a systematic change in and 
-Fdisc, similarly a change in Fi, with no corresponding change 
in -Ftot,tar would result in a reciprocal change in Fdisc. We 
therefore model the error on the term + Fdisc Vdisc in a 
Monte Carlo manner. The scaling of the Kurucz model pro- 
flle used to model the photosphere is determined by a 
minimisation over a range of scalings to flnd a best flt to 
the B and V band magnitudes of the stars as listed in the 
Hipparcos catalogue and the J H and K band magnitudes 
from the 2MASS catalogue. The percentage points of the 
distribution are used to determine a Icr limit on this scaling. 
The scaling used in the Monte Carlo modelling is taken from 
a Normal distribution with mean given by the best fitting 
scaling and standard deviation given by the level of the lo 
error. The value of Ftot is also varied between the lo error 
limits on the IRS spectrum (or rather the interpolated spec- 
trum, see above). The value of Fdisc is taken to be Ftot — F*. 
The error term on F* -I- Fdisc Vdisc is then taken from the 
standard deviation of 1000 random samplings, and is found 
to be typically at the 4% level for HD 113766 and 6% for 
HD172555 averaged over the whole wavelength range. The 
final error on the model visibility is then 
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To calculate the error on the visibility we can see from 
equation [2] that we need to add the error on the total flux, 
Ftot, tar. As we have adopted the IRS spectra for our science 
targets, this is the error on the IRS spectra which itself is 
composed of the statistical uncertainty and calibration un- 
certainty on these observations. As the spectral sampling is 
different between IRS and MIDI, we use linear interpolation 
to determine the errors at the central wavelengths of the 
spectral bins of the MIDI data. Then the flnal error on the 
visibility is given by 
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with typical values of 6% (7%) for (5Fcorr,tar/Fcorr,tar and 4% 
(3%) for 5Fcorr,tar/Fcorr,tar for HD113766 (HD172555). 

For the model visibilities equation [3] holds. We assume 
that the star has a flxed visibility of 1 (V*) and has no error. 
We also assume that Vdisc is perfectly known for any partic- 
ular model. The error on Ftot tar is discussed in the above 
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To test how well each model reproduces the observed 
visibilities we use the goodness-of-fit test. Typically this 
takes the form ^^^^(D — M)^/(7^ where A'^ is the number of 
data points, D is some observed data with associated error 
a and M is a model with no error. As both our observational 
data Vobs and model Vmod have associated errors we use the 
ratio of the observed to model visibility and compare this to 
1 (as if the model is a good flt to the data then Vmod ~ Vohs)- 
The errors for both Vohs and Vmod must be included in a. 
Our calculation then becomes 



2 

X = 



EE 



Vobs(w, n)/Vmod(w, n) 
(j(w, n) 



(15) 



Here N is the number of observations we have for a source 
and W is the number of spectral channels. To calculate a re- 
duced we divide the above by the number of data points 
(observations x spectral channels) minus the number of free 
parameters in the model (1, the FWHM, for a circularly 
symmetric Gaussian; 4, the radius, thickness, inclination 
and position angle, for a ring). If we consider equations [2] 
and |3] we can see that in dividing Vobs by Vnod we cancel 
the term in Ftot.tar. The errors from the IRS spectroscopy 
do not therefore appear explicitly in the calculation of a. 
These errors do appear implicitly in the Monte Carlo calcu- 
lation of the error on F* -I- FdiacVdisc through the calibration 
and statistical errors. The final value of a is given by the 
addition in quadrature of all the remaining error terms. As 
there is evidence of a bias in the visibility of low flux sources 
at wavelengths >12/im (see Figure [TT] and section 3.1) we 
exclude the >12/im data from the model fitting. 
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4.3 Best fitting models for HD113766 

The observed visibility functions from the MIDI observa- 
tions of HD 113766 are consistent with a partial resolution 
of the excess emission. As a first attempt to model the distri- 
bution of the excess emission we tested circularly symmet- 
ric Gaussian models with a range of FWHM 0. We tested 
1^0^ lOOmas in 60 logarithmically-spaced steps. The vis- 
ibility of the model was calculated according to the method 
described in Section 4.1, with the errors on the comparison 
of the model to the data and the goodness of fit of the model 
calculated as described in Section 4.2. The best fitting model 
is that with O=10mas, which taken over the 8-12/im range 
and the 6 baselines considered (excluding 79-F, see Section 
3.2) has a reduced of 0.32. 

Adopting an elliptical Gaussian model gives an addi- 
tional two parameters to fit. As well as the FWHM Q, 
we also consider the ellipticity of the Gaussian /, where 
vi = ucos(/), and the rotation of the model 6 which repre- 
sents the angle of the major axis of the elliptical Gaussian in 
degrees East of North. The ranges tested for each of these pa- 
rameters were Q G (1, 150)mas in 21 logarithmically-spaced 
steps, I G (0°,9Q°) in steps of 18°, and 9 G (Q°,180°) in 
steps of 10° . The best fitting model ellipticity is / = 0° , so a 
circular Gaussian model (with FWHM O = lOmas) provides 
the best fit. 

Finally, we tested ring-like distributions for the excess. 
The parameters for these models were the radius of the cen- 
ter of the ring po and the width of the ring dpo. As the Gaus- 
sian modelling showed no evidence for inclined structure, 
we do not include inclination and position angle in our pa- 
rameters for ring models. The parameter ranges tested were 
Po G (1, 100)mas and dpo G (0.2, 2.0)po (where dpo = 2po 
describes a ring from 0-2po). The best fitting model param- 
eters are po = 6mas (so ring diameter is 12mas, similar to 
the best fitting Gaussian FWHM) and dpo — 0.6po. The 
reduced for this model is 0.31. 

4.3.1 Preferred model for HD113766 

Despite the increased complexity, the ring models do not 
offer a much better fit to the data than the Gaussian mod- 
els. A simple circularly symmetric Gaussian model with a 
FWHM O — lOmas offers a good fit to the data. This is in 
keeping with the initial first order approximation to the size 
of the emitting region shown in Figure 1131 We see no evi- 
dence for a change in visibility with position angle (once the 
difference in baseline length is taken into account), and so a 
circularly symmetric distribution is to be expected to offer 
a good fit to the data. Using the percentage points of the 
distribution the 3a" limits suggest a good fit can be achieved 
with 9 < 9 < 12mas (see Figure [T6)) . The visibility of the 
best fitting Gaussian model, and the visibilities of models at 
the 5(T limits {Jd = 7mas and 14mas), are plotted on Figure 
1131 (left). We also show for comparison the visibility of the 
best fitting ring model. 

It is worth noting that for the best fitting models we 

^ If wc had instead used the VISIR spectrum as our total flux in 
the analysis of the MIDI data, the visibility functions would be 
higher (see footnote [2] The resulting best fitting Gaussian model 
would have © = 7mas. 
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Figure 16. The goodness-of-fit for circularly symmetric Gaussian 
models for the distribution of excess emission around HD 113766. 
The percentage points of the distribution are shown and indi- 
cate that only a narrow range of FWHM, 9 < < 12mas, provide 
a fit to the data within Scr. 
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Figure 17. A comparison of the observed and best fitting model 
visibilities for HD113766. A value of 1 (solid black line) would 
represent a perfect fit of the model to the data. Line colours 
and styles are as described in Figure 1121 For all baselines (ex- 
cluding 79-F) the model of a circularly symmetric Gaussian with 
FWHM = lOmas provides a good fit to the observed data (er- 
rors on Vniod/K)bs S'i'G typically around 10% across all baselines 
and the full wavelength range and are calculated as outlined in 
Section 4.2). 

achieve a reduced x^ of < 1. This would normally be termed 
an over- fit, and could be taken as evidence that the errors 
assigned to the data are too large. In this case the sources 
of error have been carefully considered in turn (see section 
4.2). A more likely possibility is that the number of degrees 
of freedom is too high. We have used the number of wave- 
length bins multiplied by the number of observations as the 
number of data-points, and subtracted the number of free 
parameters for each model considered to determine the de- 
grees of freedom. This calculation involves the implicit as- 
sumption that the data-points are independent, which for 
neighbouring wavelength bins for a particular observation 
will not be the case. This would mean that the values of 
reduced x^ are underestimated. Calcul ation of degrees-of- 
freedom in such cases is complex (see lAndrae et al.l |2010| 
for a recent discussion). A better idea of the true abso- 
lute goodness-of-fit might be achieved by examination of 
Vobs/Vmod for the best fitting model. This function is shown 
in Figure [TTl For all baselines (excluding 79-F) the circular 
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Gaussian with FWHM O = lOmas provides a good fit to tiie 
observed visibility function. 

4.4 Best fit models for HD172555 

The observed visibility functions from MIDI observations of 
IID172555 suggest a completely resolved disc on all base- 
lines. We first attempt to model the distribution of the ex- 
cess emission using circularly symmetric Gaussian models 
with a range of FWHM Q. We tested 1 sC 8 sC 200mas (the 
predicted radius of the dust is larger for this target than for 
HD113766, see Table [T| in 60 logarithmically-spaced steps. 
The visibility of the model was calculated using the method 
described in Section 4.1, with the errors on the compari- 
son of the model to the data and the goodness of fit of the 
model calculated as described in Section 4.2. The best fit- 
ting model is that with Q = 30mas, which taken over the 
8-12/im wavelength range (excluding the ozone absorption 
region) and 7 baselines (excluding problematic observations 
78-D and 83-C3, see Section 3.3) has a reduced of 1-0. 

Testing elliptical Gaussian models suggests that an el- 
liptical model provides a better fit. The best fitting parame- 
ters are found to be 6 = 48mas, I = 75° and 6 = 120°EoN. 
Such a configuration gives a slightly higher disc visibility 
along baselines at position angles close to 30°EoN, and 
slightly lower visibility at baselines close to 120°EoN. This 
orientation can be compared with Figure 1151 which shows 
that observations at baselines similar to 30° have a slightly 
higher visibility. Confidence limits on these best fitting pa- 
rameters are discussed further in subsection 4.4.1. The re- 
duced for this best fitting model is 0.36. 

We also tested ring-like distributions for the excess. The 
parameters for these models were the radius of the center of 
the ring po, the width of the ring dpo, the inclination of the 
ring to the line of sight / and the position angle of the major 
axis 6. The parameter ranges tested were po £ (6, 198)mas, 
dpo G (0.2, 2.Q)po, I G (0°, 90°) and 6* G (0°, 180°) EoN. The 
best fitting model parameters are po = 80mas, dpo = 0.5po, 
/ = 75° and 9 = 120°EoN. The reduced for this model is 
0.38. This visibility model produces higher Vinod on baselines 
close to 30° EoN and low Knod (Knod ~ Kcs) on baselines 
close to 120° EoN (compare with Figure [T^ . 

The model of [Lisse et al.. (2009.) suggests that the emis- 
sion from HD172555 might be the result of a recent massive 
collision in the system. If so, the emitting material might be 
expected to be concentrated in a clump rather than spread 
in a ring-like distribution. Dust emitting from a circumplan- 
etary region would also look like a clump orbiting the star. A 
point-like clump at an offset from the point-like star would 
have a visibility of 1 on baselines along a position angle 
90° from the direction of the offset of the clump from the 
star (e.g. if a point-like clump lay directly north of the star, 
then observations on baselines lying East- West, or at PAs of 
90°EoN, would have visibilities of 1). In the observed data 
all visibilities are <1. Along baselines in the direction of a 
point-like offset clump the visibility would be seen to os- 
cillate in a sinusoidal manner, with maxima of 1, minima 
of ~ 0.05 (averaged over the MIDI wavelength range when 
cos(27r(M(Qi — 02) + — P2))) = —1, see equation [TT|) . 
and period dependent on the distance between the clump 
and the central star. As all visibilities are compatible with 
a resolved disc at the 2a level (or better) , and no visibilities 
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Figure 18. A comparison of the observed and best fitting model 
visibilities for HD172555. A value of 1 (solid black line) would 
represent a perfect fit of the model to the data. Line colours 
and styles are as described in Figure 1141 A Gaussian of FWHM 
=30mas provides reasonable fit to the observed data according to 
the analysis (errors on Vmod/Kjbs a-re typically around 15% 
across all baselines and the full wavelength range and are cal- 
culated as outlined in Section 4.2), however a better fit can be 
found by allowing elliptical or inclined models (best fit elliptical 
Gaussian with = 48mas, /=75° and 9 = 120° EoN shown). 
An offset Gaussian clump offers no improvement over the circu- 
larly symmetric Gaussian model centered on the star. See text for 
details. 

are compatible with an unresolved source (at the > llcr level 
averaging over the 8-12pm MIDI wavelength range exclud- 
ing the 9.2-lOpm ozone absorption range), we can rule out a 
point-like offset clump as the source of the excess emission. 

To determine if a more extended offset distribution pro- 
vides a good fit to the observed visibilities, we also try to 
fit the data with a Gaussian clump offset from the central 
star. The clump is modelled in the same way as a central 
Gaussian clump, with FWHM 2 < 6 < 40mas. The clump 
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is offset from the central source (the star is again modeUed 
as a point source, so with uniform visibility of 1 everywhere) 
by a distance d, with the tested range 5 < d < SOOmas, and 
a position angle < ^ < 180°EoN. Note that we do not test 
the range 180-360° EoN as the symmetry of the visibility 
function means that the visibilities are the same for clumps 
at ^"EoN and 9 + 180°EoN. 

The best fitting clump model is for a clump of FWHM 
Q = 22mas, offset from the star by a distance d — lOmas 
at a position angle of 100° EoN. The size and offset of this 
best fitting clump model are comparable to the best fitting 
circularly symmetric Gaussian centered on the star (FWHM 
= 30mas). The offset of 100° EoN would result in higher vis- 
ibilities on baselines around ~ 10°EoN, which is what is 
seen in the observations on baseline 80-D2 (see Figure [15]). 
The reduced for this best fitting model is 1.70. Although 
this seems at first sight like a reasonable fit, as explained 
for HD 113766 the absolute values of reduced are likely to 
overstate the goodness-of-fit of the models tested here due to 
the over-estimation of the degrees of freedom in the model. 
In addition to the poorer fit, the best-fitting Gaussian clump 
model is a clump which encircles the star, and therefore is 
not a geometry which we could expect either from a circum- 
planetary di sc or from the rec ent massive collision scenario 
suggested in iLisse et all (|2009h . 



4.4.1 Preferred model for HDl 72555 

As with HD113766, the best fitting models have reduced 
< 1, which is a refiection of the over-estimation of the 
number of degrees-of-freedom. In spite of this limitation to 
our interpretation of the absolute values of reduced , we 
can see from the relative values of reduced x? that the ellipti- 
cal models and inclined rings provide a better fit to the data 
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Figure 19. The constraints we can place on the best 
fitting model parameters for an elliptical Gaussian 
model fitted to the observed visibilities of HD172555. 
We tested the goodness-of-fit for models with varying 
ellipticity (/, top left), position angle (S, top right) and 
FWHM (0 bottom left) using a Bayesian approach. 
That is, for each parameter to be fitted, we computed 
an overall value by summing the individual val- 
ues over all possible values of the other two parame- 
ters. In the bottom left we can see that any size greater 
than 35mas is favoured, which provides a lower limit 
to the size of the disc. See text for detailed discussion. 



than a circularly symmetric Gaussian or clump model. This 
is made clearer by consideration of the function V^bs/Knod 
for the best fitting models. These plots are shown in Fig- 
ure [181 It is clear from this plot that the elliptical Gaussian 
model offers a better fit than the circularly symmetric or off- 
set clump models. We find there is no significant difference 
in the goodness-of-fit offered by elliptical Gaussian and in- 
clined ring-like models. The circularly symmetric Gaussian 
and offset clump model offer similarly poor fits to the data 
from this analysis. We would expect a model with increased 
parameters to offer an improved fit to the data, and the 
fact that the offset clump model does not is reflected in the 
higher value of reduced . 

As we have found that an elliptical or inclined disc 
structure offers a better fit to the data than a circularly 
symmetric model, we want to test the limits we can place 
on the inclination/ellipticity and position angle of the disc 
with the data. We use the elliptical Gaussian model as the 
simpler model to test the parameter ranges. To determine 
the limits on the inclination of the disc model, we sum the 
values over all values of the FWHM and position angle 
for each value of the inclination I in turn. The resulting 
probability distribution is shown in Figure [191 Models with 
inclinations of 7 < 8° and I > 77° are ruled out at the 
3(T level according to the percentage points of the dis- 
tribution by this calculation, suggesting that a moderately 
inclined or elliptical disc is needed to fit the visibility data. 
We use the same technique to determine the limits on the 
other model parameters (results shown in Figure [T9]|. These 
tests show that discs with FWHM Q> < 35mas and lying 
at position angles 9 < 27°EoN or 61 > 135°EoN are also 
ruled out at the 3a level. This lower limit on the disc size 
translates to ~1AU at 29.9pc. 



© 2011 RAS, MNRAS OOP. [111221 



MIDI observations of terrestrial planet formation 21 



5 DISCUSSION 

The results presented in this paper provide the most ac- 
curate direct measurement of the dust at ~1-8AU around 
any main-sequence star. The MIDI data on HD 113766 reveal 
a partially resolved disc. The best fitting models suggest a 
disc diameter of 9-12mas (1.2-1.6AU). A circularly symmet- 
ric Gaussian provides a good fit to the data which suggests 
we do not have significant evidence for asymmetries in the 
disc structure which might be expected from a recent mas- 
sive collision for example. We also find no improvement in 
the fit if we consider elliptical Gaussians, and thus there is 
no evidence that the disc is inclined or that the binary is 
forcing an ellipticity on the disc. The binary has a spectral 
type of F5 and lies at a projected distance of 157AU. At this 
distance the binary would not be expected to have a strong 
infiuence on dust at ~1.6AU (u nless the binary has an ex - 
tremely elliptical orbit, et > O.g. lHolman fc Wiegertlll999l ). 
There is no evidence that this system has additional cold 
dust components. The Q band imaging data (which is more 
sensitive to cooler dust populations than the N band) puts a 
limit of <500mJy on cooler discs around HD113766 (see Fig- 
ure O. This hmit applies to discs with radu >0'.'15 (20AU) 
and assume a favourable disc geometry (inclined discs or 
if face-on, narrow disc geometries offer the best chance of 
detection) . 

For HD 172555 we find visibility functions that are com- 
patible with completely resolved excess emission. Modelling 
the observed data with ring-like distributions of dust sug- 
gests that the emitting material lies at a radial offset from 
the star of >35mas (>1AU at a parallax distance of 29.2pc). 
We also find tentative evidence for an inclined or elliptical 
disc structure, with major axis orientated at ~ 120°EoN. 
The TReCS imaging data on this target places an upper 
limit of 7.9AU on the extent of the disc in the N band, giv- 
ing a combined constraint of 1-7. 9AU on the dust's location. 
The TReGS Q-band imaging also indicates the presence of a 
more extended component, located at ~8AU from the star. 
This component is inclined or elliptical in structure, lying at 
a position angle of ~110°EoN. The N band imaging shows 
that the 10/^m emission must come from a region inside this, 
and from the MIDI results we infer a similar postion angle to 
that seen in the Q band image, but at a low significance. This 
suggests that the two may be linked in an extended disc, 
with the N and Q band data probing the inner and outer 
regions of the distribution respectively. However, HD 172555 
may also have two spatially distinct com ponents, similar to 



other multi-component d iscs (e.g. - q Tel. iSmith et al] 2009al : 
/3 Leo. IStock et al.ll20ld : ?? Corvi. ISmith et al.ll2009"bl ). HD 



17255 5 could be viewed a s another young solar system (see 
rj Tel : ISmith et al.ll2009al ). and so could alternatively be in- 
terpreted as having two debris belts somewhat analogous to 
the asteroid and Kuiper belts in the Solar System. Neither 
configuratio n rules out the rec ent massive collision scenario 
suggested bv lLisse et al] (| 2009h for the origin of the emission. 
An alternative extended disc distribution could also explain 
the results presented in this paper, with the N and Q band 
data probing the inner and outer regions of the distribution 
respectively. Again, such a distribution does not rule out 
the possibility of a recent massive collision, although if the 
majority of the emission is believed to have arisen from a 
collision, dynamical models will be required to confirm that 



the resulting ejecta can be spread over such a wide spatial 
range. 

The young ages of the systems studied here (16Myr 
and 12Myr), combined with the lack of observed cold dust 
in these systems and the location from which the observed 
emission is seen to arise all sugges t we are wi tnessing on- 
going terrestrial planet formation. iLisse et al.l (f2009 ) con- 
cluded that the spectrum of HD 113766 is consistent with 
emission from a disrupted S-type asteroid (of ^ 320km in 
radius) or an asteroid belt made up of a large number of 
small S-type asteroids. Such a belt is a natural consequence 
of ter restrial planet f ormation (s e e e.g. , lKenvon fc BromlevI 
l2004h . For HD172555 rLisse et al.l (|2009l ') suggest that the ex- 
cess emission arises from a recent massive collision similar to 
that believed to have formed the Earth-Moon system. If this 
is a true model of the source of the excess then we might ex- 
pect a clumpy dust distribution, however we do not see any 
evidence for this in the visibility functions. IWvatt fc DentI 
('2002') presented a model for the evolution of a clump follow- 
ing a massive collision in a debris disc, which they found to 
be an unlikely source for the asymmetric dust distribution 
in Fomalhaut. Assuming a radius of 5.8AU, a parent plan- 
etesimal radi us of 1000km wit h density 2.5 g cm~^ (from 
the fitting of Lisse et al. 200!^) and a specific incident en- 
ergy required for catastrophic destruction of a planetesimal 
(5J)=200 J kg~^, we find that the time taken for a dust 
cloud produced in a massive collision to occupy half the disc 
(and thus no longer be observable as a clump) is ~200 years. 
As the constraints from the modelling of the IRS spectrum 
on the lifetime of the dust arising from a massive collision 
suggest that the collision should have occurred in the last 
O.lMyr, clump- like emission would only be expected to have 
lasted for the first 0.2% of the lifetime of the emission (this 
is consistent with a clum p lifetime of 20-30 orbital periods 
seen in the asteroid belt. iMichel et a"l]|200ll ). It is not sur- 
prising therefore that the emission we observe is not seen 
to be clump-like if the origin of the emission is a massive 
collision. The poor fit to the observed visibilities provided 
by clump models suggests that we can rule out the possi- 
bility of the emission arising from circumplanetary regions 
(around an as yet undetected planet). 

The constraints we can place on the radial location of 
the emitting dust around HD113766 and HD172555 with 
these observations are of particular interest, as when com- 
bined with spectral analysis they allow the possibility of 
constraining models for the dust grains themselves. This is 
important because predictions of the location of the debris 
based on spectral analyses alone can be highly dependent 
on assumptions made about the grain properties such as 
size. For HD113766, the location of the du st inferred from 
spectral analysis is 13mas (1.7AU lLisse et a l. 2008), slightly 
larger than our MIDI observations imply. For HD172555 
dust has been pre dicted to arise from a region around 5.8AU 
(|Lisse et aLlbOOS^ *). at the larger end of the 1-7.9 AU scales 
that the MIDI and TReCS observations presented here im- 
ply. We believe these discrepancies are due to the modelling 
methods used, scaling observations made during the Deep 
Impact Temple 1 experiment to account for the luminosity of 
each target star. As our new interferometric data constrains 
the dust in both system s to be loc a ted in regio ns smaller 
than p redicted from the lLisse et al.l (|2008l ) and iLisse et al.l 
models, the implication is that the dust composi- 



tnan p : 
(|2009l ) 
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tion is different from the dust released in the Deep Impact 
event such that it would be hotter if placed at the same dis- 
tance from any given star. These constraints will allow the 
next-generation models to treat the dust composition and 
location in a self-consistent way. 

In this paper we have presented new VISIR spec- 
troscopic and MIDI interferometric observations of the 
youngest debris disc hosts with emission on <C10AU scales, 
HD113766 and HD172555. We have additionally presented 
new VISIR imaging observations of HD113766, which rule 
out large-scale high surface brightness emission. We find that 
the visibility functions measured for HD 113766 are consis- 
tent with symmetric emission in the region between 0.6- 
0.8AU (or a diameter of 1.2-1.6AU). For HD172555 the vis- 
ibility functions observed suggest that the emission is com- 
pletely resolved with the VLTI, placing a lower limit of lAU 
on the spatial extent of the emitting regi on. A new anal- 
ysis o f TReCS data originally presented in lMoerchen et al] 
has revealed extended emission around HD 172555 in 
the Q band, consistent with an inclined disc at a median ra- 
dius of 7.9AU from the star. This means that the HD172555 
system could contain multiple disc components, or that we 
are observing different regions of a broad dust distribution 
around this stars in the different bands. These observations 
are the first that resolve the mid-infrared excess emission 
around these targets. In both cases the observational data 
are compatible with current models that suggest we are ob- 
serving two systems in the midst of terrestrial planet forming 
processes. 
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